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Abstract.  We  describe  a  new  shortest  paths  algorithm.  Our  algorithm  achieves  the 
same  0(nm)  worst-case  time  bound  as  Belbnan-Ford  cilgorithm  but  is  superior  in  prac¬ 
tice. 


1.  Introduction 

The  BeUman-Ford  algorithm  [1,  4,  7]  is  a  classical  algorithm  for  the  single-source  shortest 
paths  problem.  The  5dgorithm  runs  in  0{nm)  time  on  a  graph  with  n  nodes  and  m  arcs.  This 
is  the  best  currently  known  strongly  polynomial  bound  for  the  problem  (see  [6]  for  the  best 
weakly  polynomial  bound). 

In  practice,  however,  the  BeUman-Ford  algorithm  is  usually  outperformed  by  the  deque 
.algorithm  of  D’Escopo-Pape  [9]  and  by  the  two-queue  algorithm  of  Pallottino  [8].  The  work- 
case  time  bounds  for  these  algorithms,  however,  are  worse  than  those  for  the  BeUman-Ford 
algorithm.  The  deque  algorithm  may  take  exponential  time  in  the  worst  case  [10]  and  the 
two-queue  algorithm  may  take  time. 

We  propose  a  new  topological-scan  algorithm  for  the  shortest  paths  problem.  The  algorithm 
combines,  in  a  natural  way,  the  ideas  of  the  BeUman-Ford  algorithm  with  those  for  finding 
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procedure  SCAN(t;); 

for  all  (v,iy)  E:  E  do 

ifd{v)  4-  i{vj  u;)  <  d(iu)then 
d{w)  4—  d{v)  H-  £(v,  w); 
S{w)  ^  labeled; 

7r(u;)  ^  V] 

5(t;)  scanned; 
end. 


Figure  1.  The  scan  operation. 

shortest  paths  in  acyclic  graphs  (see  e.g.  [2]).  The  algorithm  achieves  the  same  O(nTn)  time 
bound  as  the  Bellman- Ford  algorithm,  but  outperforms  it  on  most  problem  instances.  A 
preliminary  study  shows  that  the  cdgorithm  is  competitive  with  the  deque  and  the  two-queue 
algorithms. 


2.  Definitions  and  Notation 

The  input  to  the  single-source  shortest  paths  problem  is  (G,s,/),  where  G  =  {V,E)  is  a 
directed  graph,  /  :  R  is  a  length  function,  and  s  G  V  is  the  source  node.  The  goal  is  to 

find  shortest  paths  from  s  to  all  other  nodes  of  G  or  to  find  a  negative  length  cycle  in  G.  If 
G  has  a  negative  length  cycle,  we  say  that  the  problem  is  infeasible.  We  assume,  without  loss 
of  generality,  that  all  nodes  are  reachable  from  s  in  G  and  that  G  has  no  multiple  arcs.  The 
latter  assumption  edlows  us  to  refer  to  an  arc  by  its  endpoints  without  ambiguity.  We  denote 
|V|  by  n  Eind  \E\  by  m. 

A  potential  function  is  a  real- valued  function  on  nodes.  Given  a  potential  function  d,  we 
define  a  reduced  cost  function  :  E  K  hy 

id{v,  w)  =  £{v,w)  +  d(v)  —  d(w). 

Given  a  potential  function  d,  we  say  that  an  arc  a  is  admissible  if  /^(a)  <  0,  and  denote  the 
set  of  admissible  arcs  by  E^.  The  admissible  graph  is  defined  by  Gd  =  {V,Ed). 

A  shortest  paths  tree  of  G  is  a  spanniug  tree  rooted  at  s  such  that  for  any  «  G  V,  the  reversal 
of  the  v  to  s  path  in  the  tree  is  a  shortest  path  from  s  to  v. 

3.  The  Label- Correcting  Method 

Our  edgorithm  is  based  on  the  general  label-correcting  method  for  solving  shortest  paths 
problem.  In  this  section  we  briefly  outline  the  method.  (See  e.g.  [2,  11]  for  more  detail.) 

For  every  node  v,  the  method  maintains  its  potential  d[v),  pment  7r(u),  and  status  S{v)  G 
{unreached,  labeled,  scanned}.  Initially  for  every  node  v,  p(v)  =  oo,  ir(t;)  =  nil,  and  5(v)  <- 
unreached.  The  method  starts  by  setting  d(s)  *-  0  and  S(s)  =  labeled,  and  applies  the  scan 
operation  to  labeled  nodes  imtil  none  exists,  in  which  case  the  method  terminates. 
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The  SCAN  operation  applies  to  a  labeled  node  v.  The  operation  is  described  in  Figure  1. 
Note  that  if  v  is  labeled,  then  d{v)  <  oo  zind  d{v)  +  l{v,w)  is  finite. 

The  method  terminates  if  atnd  only  if  G  does  not  have  negative  length  cycles.  If  the  method 
terminates,  the  parent  pointers  define  a  correct  shortest  paths  tree  and,  for  amy  v  eV,  d{v)  is 
the  shortest  path  distemce  from  s  to  t;.  The  label- correcting  method  can  be  easily  modified  so 
that  if  G  has  negative  cycles,  the  method  finds  such  a  cycle  and  terminates. 

Different  strategies  for  selecting  labeled  nodes  to  be  scanned  lead  to  different  algorithms. 
The  BeUman-Ford  algorithm  medntains  the  set  of  labeled  nodes  in  a  FIFO  queue.  The  next 
node  to  be  scanned  is  removed  from  the  head  of  the  queue,  a  node  that  becomes  labeled  is 
added  to  the  tsul  of  the  queue. 

The  D’Escopo-Pape  algorithm  uses  a  deque  (a  queue  that  edlows  insertions  at  either  end) 
to  m£iintain  the  labeled  nodes.  The  next  node  to  be  scEumed  is  removed  from  the  head  of  the 
queue.  A  node  that  becomes  labeled  is  added  to  the  tail  of  the  queue  if  this  is  the  first  time 
this  node  became  labeled,  and  to  the  head  of  the  queue  otherwise  (t.e.,  if  it  has  been  scEinned 
previously). 

Pallottino’s  two-queue  algorithm  uses  two  FIFO  queues,  one  of  low  priority  and  another  of 
high  priority.  Each  labeled  node  appears  on  exactly  one  of  the  queues.  The  next  node  to  be 
scaimed  is  removed  from  the  head  of  the  high  priority  queue  if  the  queue  is  not  empty  and 
from  the  head  of  the  low  priority  queue  otherwise.  A  node  that  becomes  labeled  is  added  to 
the  tail  of  the  low  priority  queue  if  this  is  the  first  time  this  node  became  labeled,  and  to  the 
tail  of  the  high  priority  queue  otherwise. 

Dijkstra’s  Eilgorithm  [3]  maintains  the  set  of  labeled  nodes  in  a  priority  queue  eind  selects  a 
labeled  node  with  the  minimum  potential  as  the  next  node  to  be  scanned.  When  the  length 
function  is  nonnegative,  this  algorithm  nms  in  0(m  +  nlogn)  time  if  Fibonacci  heaps  [5]  are 
used  to  implement  the  priority  queue.  If  the  graph  has  negative  length  Etrcs,  the  algorithm  can 
take  exponential  time. 


4.  The  Topological-Scan  Algorithm 

Suppose  both  v  and  w  are  labeled  Emd  ti{v,w)  <  0.  Then  it  is  better  to  scan  v  before 
w,  since  we  know  that  when  we  scan  v,  d(w)  will  improve  and  w  will  become  labeled.  The 
topological-scan  algorithm  uses  a  generalization  of  this  idea. 

To  simplify  the  algorithm  description,  we  assume  that  G  has  no  cycles  of  length  zero  or  less, 
and  therefore  for  any  d,  Gd  is  acyclic.  Later  we  show  how  to  get  rid  of  this  assumption. 

The  topological-scan  algorithm  maintains  the  set  of  labeled  nodes  in  two  sets,  A  and  B. 
Each  labeled  node  is  in  exactly  one  set.  Initially  A  =  0  and  B  =  {«}.  At  the  beginning  of 
each  pass,  the  algorithm  uses  the  set  B  to  compute  the  set  A  of  nodes  to  be  scanned  during 
the  pass,  and  resets  B  to  the  empty  set.  A  is  a  linearly  ordered  set.  During  the  pEiss,  elements 
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axe  removed  according  to  the  ordering  of  A  and  scauaned.  The  newly  created  labeled  nodes  are 
added  to  5.  A  pass  ends  when  A  becomes  empty.  The  algorithm  terminates  when  B  is  empty 
at  the  end  of  a  pass. 

The  algorithm  computes  A  from  B  as  follows. 

(1)  For  every  r  £  5  that  has  no  outgoing  arc  with  negative  reduced  cost,  delete  v  from  B 
and  mark  it  as  scetnned. 

(2)  Let  A  be  the  set  of  nodes  reachable  from  B  in  Gd-  Mark  aU  nodes  in  A  as  labeled. 

(3)  Apply  topological  sort  to  order  A  so  that  for  every  pair  of  nodes  v  and  w  in  A  such 
that  (u,  w)  £  Gd,  V  precedes  w  smd  therefore  v  will  be  scanned  before  w. 

Theorem  4.1.  The  topological-scan  algorithm  is  correct 

Proof.  Note  that  if  a  node  v  has  no  outgoing  arc  with  negative  reduced  cost,  then  the  only 
side-effect  of  scan(v)  will  be  the  change  of  the  status  of  v  to  scanned.  Note  dso  that  the 
label- correcting  method  remains  correct  if  we  mzdce  a  scanned  node  into  a  labeled  node,  and 
that  all  nodes  reachable  from  5  in  Gd  at  step  2  have  finite  potentials  and  therefore  are  either 
labeled  or  scanned.  Finally  recall  that  we  assumed  that  Gd  is  acyclic,  and  thus  step  3  is  well 
defined.  These  observations  and  the  fact  that  the  label-correcting  method  is  correct  imply  that 
the  edgorithm  is  correct.  ■ 

Analysis  of  the  topological-scan  algorithm  is  very  similar  to  that  of  the  BeUmem-Ford  algo¬ 
rithm,  as  given  for  example  in  [2,  11]. 

Theorem  4.2.  The  topological-scan  algorithm  runs  in  0{nm)  time. 

Now  suppose  G  has  cycles  of  zero  or  negative  length.  In  this  case  Gd  need  not  be  acyclic. 
If,  however,  Gd  has  a  negative  length  cycle,  we  can  terminate  the  computation.  If  Gd  has  zero 
length  cycles,  we  can  contract  such  cycles  eind  continue  the  computation.  This  can  be  easily 
done  while  maintaining  the  0{nm)  time  bound.  (See  e.g.  [6].) 

5.  Concluding  Remarks 

To  obtain  a  practical  implementation  of  the  algorithm,  it  is  important  to  implement  the 
procedure  of  constructing  A  from  B  efficiently.  This  construction  can  be  done  in  one  depth- 
first  search  computation.  It  uses  stacks  to  implement  A  and  B.  The  computation  popes  nodes 
from  B  one  by  one.  For  each  removed  node  v,  if  w  was  already  visited  by  the  current  depth-first 
search,  we  do  nothing.  If  u  was  not  visited  and  v  has  no  outgoing  arcs  of  negative  reduced  cost, 
we  mark  v  as  scanned.  If  v  was  not  visited  and  has  an  outgoing  arc  of  negative  reduced  cost, 
we  visit  in  the  depth-first  order  aU  nodes  reachable  from  u  in  Gd  md  not  visited  previously. 
At  the  end  of  the  depth-first  search  visit  of  a  node,  the  node  is  marked  labeled  and  pushed 
on  A.  Note  that  this  procedure  can  be  implemented  so  that  it  looks  only  at  the  nodes  which 
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were  on  B  at  the  beginning  of  the  computation  or  zu-e  on  A  at  the  end  and  each  such  node  is 
exsunined  exactly  once.  Thus  if  the  average  number  of  labeled  nodes  at  the  beginning  of  every 
phase  is  small,  the  average  time  per  phase  is  small. 

An  alternative  way  of  deeding  with  potential  cycles  in  Gj  is  to  ignore  some  arcs  of  Gj  and 
produce  a  topological  ordering  with  respect  to  the  (acyclic)  graph  G'  obtained  from  Gd  by 
deleting  the  ignored  arcs.  A  stemdard  way  to  implement  topologiced  sort  is  by  using  depth-first 
seeirch.  If  we  ignore  the  arcs  that  depth-first  search  classifies  as  back  arcs  (see  e.g.  [2]),  then  we 
obtadn  a  topological  ordering  of  nodes  with  respect  to  G'.  To  assure  that  the  algorithm  detects 
negative  cycles,  we  can  use  standard  techniques  developed  for  the  label-correcting  method. 
The  resulting  algorithm  still  runs  in  0{nm)  time. 

Finally,  we  would  like  to  note  that  the  topological-sczm  algorithm  runs  in  0(7n)  time  if  the 
input  graph  G  is  acyclic. 
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